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Abstract. We performed detailed investigation into cosmological perturbation of f(T) theory 
of gravity coupled with scalar field. Our work emphasizes on the exact way of gauge fixing and 
we examined all possible modes of perturbations up to second order. This includes in addition 
to the usual scalar, vector, and tensor modes, also pseudoscalar and pseudovector modes; 
although we find that there is no gravitational propagating degrees of freedom in the scalar, 
pseudoscalar, vector, as well as pseudovector modes. We also find that the scalar and tensor 
perturbations have exactly the same form as their counterparts in usual general relativity 
with scalar field, except that the factor of reduced Planck mass squared M^^ = 1/{8ttG) 
that occurs in the latter has now been replaced by an effective time- dependent gravitational 
coupling —2{df / dT)\T=To : with Tq being the background torsion scalar. The absence of 
extra degrees of freedom of /(T) gravity at second order linear perturbation indicates that 
/(T) gravity is highly nonlinear, and one cannot conclusively analyze stability of the theory 
without performing nonlinear analysis that can reveal the propagation of the extra degrees 
of freedom. 
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1 Introduction 

Einstein proposed the idea of teleparallelism, or Fernparallelismus, to unify gravity and elec- 
tromagnetism into a unified field theory in 1928 [1]. Unlike general relativity (GR) in which 
the Levi-Civita connection gives rise to curvature but vanishing torsion, in teleparallelism, 
spacetime is endowed with a connection with vanishing curvature, but nonzero torsion. Since 
the curvature is identically zero, parallel transport of a vector is independent of path. This 
is the origin of the name teZeparallel, which means "parallel at a distance" (Einstein's quest 
for unified field theory via teleparallelism is an interesting history and can be found in [2]). 
It has since been established that general relativity can in fact be re-casted into teleparallel 
language [3] [4] [5] [6] , known as the Teleparallel Equivalent of General Relativity (TEGR) . For 
an interesting formulation of TEGR as higher gauge theory, see [7]. Due to the need to 
understand the acceleration of the universe, various theories of modified gravity had been 
introduced, among which is the attempt to generalize TEGR to f(T) theory of modified 
gravity in the same spirit as generalizing general relativity to f{R) gravity [8] [9]. 

Given the FLRW (Friedmann-Lemaitre- Robertson- Walker) solution with suitable choice 
of the scale factor that describes an accelerating universe, the next step is to understand the 
evolution of structures on the FLRW background. From the observational point of view, this 
has been discussed in [10]. On the theoretical level, perturbation theory is a useful tool for 
this purpose as it can reveal some properties of the dynamical modes of a gravitational theory. 
In the case of f{T) theory, Li et. al [11] had shown that there are generically, 5 degrees of 
freedom in f{T) gravity. Comparing with GR, which has only 2 degrees of freedom, there are 
thus 3 extra degrees of freedom, which the authors suggest could correspond to either one 
massive vector field or one massless vector field with one scalar field. Note that due to the 
high symmetry of FLRW metric, it can be checked that there is no extra degree of freedom 
at the background level. This is actually to be expected: When restricted to some particular 
symmetry such as FLRW geometry or spherical geometry, some degrees of freedom must 
vanish by symmetry. In particular spatial isotropy would not allow a non-vanishing spatial 
vector field, which selects a particular special direction. At the level of perturbation however, 
one must carefully check if any extra degree of freedom manifests. In the case of f{T) gravity. 
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as pointed out in [11], a number of extra degrees of freedom do not appear in Minkowski 
spacetime, even at linear perturbation level. 

We remark that extra propagational degrees of freedom can sometimes create additional 
interaction between two bodies. In order for the theory to be consistent with observational 
data, we need some mechanisms such as the suppression of such interaction or nonlinear 
effect that prevents the propagations of the extra modes. If f{T) gravity is the exact theory 
describing our universe, then FLRW solution should be stable. However stability can only 
be established after we know the behaviors of all the degrees of freedom of f{T). 

In order to address these questions, we analyze the perturbative properties on FLRW 
background in f{T) gravity. Although linear perturbation on this particular background had 
been previously investigated in various papers [12] [13] [14], in some of these a few equations 
are ignored, and in others the way of carrying out the decomposition is incomplete. We 
believe that it is thus worth while to re-visit the problem by considering more general type 
of perturbation, including all possible modes: scalar, pseudoscalar, vector, pseudovector, 
and tensor, by emphasizing on the exact way of gauge fixing. Furthermore, we perform our 
analysis up to second order of action in perturbation, which is important for checking the 
ghost condition. 

The rest of this paper is organized as follows. For the sake of self-completeness, in 
Sec. 2 we shall review Teleparallel Equivalent of General Relativity and f(T) theory. Sec. 3 
is our analysis: In Subsec.3.1, we show the exact way of the decomposition, followed by the 
analysis of the second order action in each mode from Subsec.3.2 to 3.5. Sec. 4 is devoted to 
the summary of this paper and discussion. 

2 Prom Teleparallel Equivalent of General Relativity to /(T): A Review 

As in general relativity, the underlying structure of TEGR is a spacetime manifold M. At 
every point p £ M in local coordinate chart {x^, • • • ,2;"}, the tangent space TpM at p is 
spanned by the coordinate vector fields {di,--- ,dn}- The set jdx^,--- forms the 

basis for the dual space T*M. Specializing to 4-dimension, the tangent space is equipped 
with Lorentzian metric g with signature (— , -|-, + , -|-). We will label all spacetime coordinates 
by Greek indices that run from to 3, with denoting the time dimension, while all spatial 
coordinates will be labeled by i,j, k, . . . that run from 1 to 3. 

Let {e^ix)} be an arbitrary basis of TpM. We can express the total covariant derivative 
VeA as: 

VeA{x) = r^Ai^)eB{x), (2.1) 
where the coefficient is the connection 1-form, satisfying 

r^^(x) = (^?^, Ve^) = r%{x)dx'^, (2.2) 

where |t?^(x)} is the dual frame or coframe of {eB{x)} and {dx^^} is the dual basis of 
coordinate basis {d^}. While Greek indices correspond to coordinate basis, Latin indices 
refer to the so-called "internal" indices, which correspond to non-coordinate frame in the 
fiber. That is, eA = ^^a^i^ '^^ ~ ^\dx^. In this work, A, C, ... take the values 0, 1, 2, 3 
while a, 6, c, ... take the value 1, 2, 3. 

Assuming that spacetime is parallelizable (i.e. there exist n vector fields {vi, ...,Vn\ 
such that at any point p £ M the tangent vectors flip's provide a basis of the tangent 
space at p), we can view the mapping between the bases in coordinate frame {d^} to that of 
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non-coordinate frame {e^} as an isomorphism TM — )• M x M^. This follows from the fact 
that any parallelizable n-dimensional manifolds have trivializable tangent bundle (i.e. can 
be written as the direct product of M and M"). Note that the frame field depends only on 
the affine structure of the manifold and hence a priori has nothing to do with the metric. 
We now relate the frames with the metric by equipping the space M"* with a flat Minkowski 
metric r/^^ such that 

9^lu = ??Ase^M^^^'■ (2.3) 
In other words, on the manifold M, we choose arbitrarily a frame ca at each point, locally 
on some open chart U C M (by parallelizability, this can be extended globally), and define 
a metric 77 on [/ by 

r]{eA,eB) =Vab- (2-4) 

This imposes orthonormality on the tetrads. 

We remark that for gravity in (3+l)-dimension, things work out nicely because of Steen- 
rod's theorem, which guarantees that all orientable 3-manifolds are parallelizable [15] [16]. 
Thus a 4-dimensional spacetime with orientable spatial section is also parallelizable. To put it 
slightly differently, if we believe that any spatial slice of spacetime is an orientable 3-manifold 
(and as such parallelizable) with some initial data that can be propagated uniquely in time, 
in the manner of 3+1 decomposition of ADM [17], then we should believe that the entire 
spacetime is parallelizable. We remark also that a non-compact 4-dimensional Lorentzian 
manifold M (all spacetime is non-compact, for otherwise causality will be violated) admits 
a spin structure if and only if it is parallelizable, a result known as Geroch's theorem [18]. 

w 

One can now introduce the Weitzenbock connection [19] [20] V defined by 

VxY = {XY^)eA, (2.5) 
where Y = Y^ca- That is, we require that the tetrads are teleparallel, i.e. covariantly 

w 

constant: Vx^A = 0. This amounts to defining the connection coefficient by 

w 

r V = ^Adue""^ = -e\d,el (2.6) 

Given any vector field X, there always exists such connection that gives rise to parallelization 
on M, assuming of course M is parallelizable. Furthermore, the connection is unique [21]. 
The torsion tensor is defined by 

W WW 

T{X,Y) = VxY -VyX -[X,Y] (2.7) 
= X^Y^[eA,eB]. (2.8) 

This expression is generically not zero since the basis vectors {e^} are not integrable in 
general. In local coordinates, the components of the torsion tensor is 

w w w 

T\u ^ r^,^ - = e^{d^e\ - d,e\) / 0. (2.9) 

The Weitzenbock connection is curvature-free. One can verify this either by direct 
(but tedious) calculation in coordinates, or simply by noting that the definition of Riemann 
curvature endomorphism 

w / w w w w w \ 

R{X, Y)Z = f VxVy - Vy Vx - V[x,y] j Z, (2.10) 
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and V x^A = {X5\)ec = together implies 

W WW WW w 

R{eA,eB)ec = Ve^(Ve3ec) - Ves(Ve^ec) - V[e^_eB]ec, (2.11) 

every term of which vanishes identicahy and so the Riemann curvature tensor is equal to 
zero. We would like to emphasize that curvature is not a property of the manifold itself, but 
depends on the choice of connection. 

The difference between the Weitzenbock connection and the Levi-Civita connection is 
the contortion tensor^ 

K\, = -iT,\ + T/,-T\A, (2.12) 

or equivalently, 

W 1 / w \ 

K''% = - 2 ( - ^''^P - ^p'"' ) ■ (2-13) 

The Einstein-Hilbert action of general relativity (here the Ricci scalar R is with respect 
to the Levi-Civita connection) 

S=^ [ d^x./^R, K = SttG, (2.14) 



2k 

can be written as 



1 

'2k 



d'^x eT, (2.15) 
where e = | det(e'^)|, which is equal to \/—g in GR, and 

www 

T = Sp^'TP^,, (2.16) 



with 



w '\ / w w w \ 

V'^ ^ - (k^% + 6>;t-\ - 5^r"^J . (2.17) 



w 



The quantity T is the so-called torsion scalar. Explicitly, 

W '] W W '\ W W WW 

T = -r%^V^ + -T^r'?^ - Tp/T''^',. (2.18) 
The equivalence between Eq.(2.14) and Eq.(2.15) comes from the equality 

w w 

i2 = -r-2V''r%, (2.19) 

where V is the covariant derivative in GR. 

The Einstein-Hilbert action in the form of Eq.(2.15) is known as the TEGR action. 
f{T) gravity is obtained by replacing T in TEGR action with a function f{T). 

Note that unlike TEGR, we do not have local Lorentz invariance in /(T) gravity, be- 

w 

cause the total derivative term VT'^^j^ in Eq.(2.19) breaks local Lorentz invariance (see e.g. 
discussions in [23] and [24] ) . Thus unlike in general relativity in which we can change coordi- 
nate systems and use any frame fields, we have no such freedom in f{T) gravity. In particular, 
this implies that we cannot directly extract the tetrad from the metric in a straightforward 



'^Note the proper term is contortion instead of contorsion. See e.g. [22]. 
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manner. TEGR is special since its action does not determine the admissible frame but only 
the metric. Nevertheless, in the case of FLRW cosmology with flat spatial section, 

ds^ = gf.^dxf'dx'' = -dt^ + a'^6ijdx'dx\ (2.20) 

where a = a{t) is the expansion factor, the standard choice of tetrads, or equivalently the 
cotetrads, related to the coordinate one-forms by 

gO^dx^ = dt, (2.21) 
e'^^dx'' = a5atdx\ (2.22) 

is a suitable choice of frame, albeit not the most general one (for example, one may wish 
to consider spacetime with less symmetry such as Bianchi models. In fact this will be quite 
important for the analysis of /(T) gravity. We will further comment on this issue in the 
Discussion section). 

From now onwards, we will suppress all explicit overscript '-u;' on the torsion scalar and 
connection coefficients etc. 



3 f{T) Gravity and Its Cosmological Perturbation 

In this work, we consider the cotetrads 

^^dx^" = dt, (3.1) 
e'^^dx'' = a5aidx\ (3.2) 

that is, the background metric is taken to be the usual FLRW geometry with flat spatial 
section given by Eq.(2.20). We will denote the Hubble parameter hy H = a/ a. 
The gravitational action is 

Sg = j d^x ef{T). (3.3) 

We also introduce scalar field into our theory as a toy model for the matter sector. The 
matter scalar field is described by the action 

S^ = j d^x e (-]^g^'^d^(t)d,(t> - y (</.)) , (3.4) 

where V{(j)) is an unspecified potential. 

3.1 Decomposition of Degrees of Freedom and Gauge Freedom 

Using the background symmetry, we can decompose the perturbative cotetrad into five 
scalars, one pseudoscalar, three vectors, one pseudovector and one tensor^. In the following 
calculations, all equality is understood to be either exact, or exact up to second order only. 



^ In [14], the authors had discussed the decomposition but it is imperfect. Although they had tried to 
write the antisymmetric part of (Je"; in terms of vector and scalar, this should more appropriately be written 
in terms of pseudoscalar and pseudovector. While pseudovector can be written in terms of vector due to the 
relation Eq.(3.9), pseudoscalar can never be written in terms of scalar. 
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The exact forms of decomposition are 

6e't = f , (3.5) 

= a {diP + Ui) , (3.6) 

Se^'t = 5a^ {d,B + vi) , (3.7) 



(5e" = a6aj 



6ijij + didjE + diWj + djWi + hij + eijk ( 5^(7 + Vu 



(3.8) 



where /3, V and E are scalars, o" is a pseudoscalar, Uj, and tfj are vectors, 1^ is a 
pseudovector and hij is a tensor. We note that vectors and pseudovectors can be coupled 
with each other in the following form: 

f-ijk (diUj) Vfc- (3.9) 

If the theory preserves parity, then different modes of decomposition cannot be coupled with 
each other, except for the coupling between vectors and pseudovectors. This means that 
we can, with the exception of vector and pseudovector modes, study the different modes 
separately. 

As in GR, the action of f{T) gravity is invariant under coordinate transformation: 

xt' ^ x"" = xf" + ^f'{x). (3.10) 
Under this transformation, the cotetrad transforms as 

Ki-') = -i^Stt- (3-11) 

This leads to the gauge transformation of perturbative variables on FLRW metric (2.20) as 

$^$' = ^.-^*. (3.12) 
(3^f3' = (3--e, (3.13) 



a 



a 



a 



B^B' = B-l^--C], (3.14) 



^^^' = ^--^*, (3.15) 
a 

E ^ E' = E - -C, (3.16) 
a 

a^a' = a, (3.17) 

Ui — n- = u'^, (3.18) 

^ v'^ = v^- (ef ) - , (3.19) 

Wi ^ w'i = Wi - (3.20) 



2a 
1 



Vi^V^ = Vi--eijkdjC'^^\ (3.21) 



a 



^ h[j = hj, (3.22) 



where we have decomposed ^* as ^* = o ^{di^ + Sjf^) and dot denotes the partial derivative 
with respect to time. We know we can choose the gauge where one of /3 and ip, one of 
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B and E and one of Vi, Wi and Vi vanish. That is, by gauge-fixing, Vi can be set to zero 
and thus we only need to consider perturbative analysis of vector modes. We now begin 
our perturbative analysis of the various modes of degrees of freedom, starting from scalar 
perturbation. 

3.2 Scalar Perturbation 

We choose the gauge where P = B = 0. For simplicity, instead of ^, we use another variable 
C, which is defined as 



C = V' + 



AE 



(3.23) 



where A = Sijdidj. Then, cotetrad and metric are written as 



dx^' = (1 + dt, 



e^dx'' = 5a^ |a 5ij{l + C) + (^d,dj - ^A6,,^ E | 



dx^ , 



.. d 



a~^iy5a^{l - C + C') - 5aj (d^Oj - E + 26ajC [did, - E 



+5, 



didk - -A5ik 1 E 



g^udx^dx'' = -{l + ^fdt^ 



+a'{ 5ij{l + 0' + 2(1 + C) didj - -AS,, E 



djdk - -A6jk 1 E 



1 



_d_ 



+ 



d.dk - -A6.k 1 E 



djdk - ^A(5jfc ) E 



dx^dx^ , 



^ d_ 

Ox'* dx'' 



[1 - 2$ + 3$^ 



d_d_ 

^dtdi 



+a-^!^5,j - 26^jC - 2 [d^dj - ^A6^j^ E + 36ijC^ + 6C (d^d, 



+3 



didk - -A5ik 1 E 



djdk - ) E 



d d 
dx^ dx^ 



(3.24) 
(3.25) 

(3.26) 
(3.27) 



(3.28) 



^A6ij 1 E 



(3.29) 



Then, we have 



e = a 



1 + $ + 3C + 3$C + 3C' + li^Ef - \{didjE){didjE) 



(3.30) 
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The components of the torsion tensors are 



tj 



T 



jk 



0, 



H5ij + 5ijC + { d,dj - -^5ij ) E - {dAE) 



dkdj - -A5kj ) E 



mdjE) + (-C + Iae) 



S^jC + { didj - -AS,, ) E 



(3.31) 
(3.32) 



(3.33) 



dj6,kC - lAd,6ikE + (^-C + ^AE^ (^dj6ikC - ^Aa,<5,fci?^ + ^{didkE){AdjE) 
-{didkE){djO - (j o k), (3.34) 



T^t^ = 3H + 3C- 3CC - {didj E) {didj E) + -{AE){AE), 



+ 2diC - lAdiE - <^d,<^> + (^-C + ^A^^ (^2diC - l^9iE^ 
+^{AE){AdiE) - ^{did,E){djAE) - {AE){d,C) + {djdiE){djC), 



(3.35) 



(3.36) 



where (j -H- k) in Eq.(3.34) denotes the permutations of the all previous terms with respect 
to j and k. 

Define = T"^^. Then, we have the following: 



TaT^gf"" = -m'^ - 18HC + 18H'^^ - 9C^ + 36H^C - 27H'^^'^ + 18HCC 



r2/F>2 



+6H {didj E) {didj E) - 2{AE){AE)H + a"' 
= -3H^ - QHC + 6H'^<^ + a-^{di^f + 2o" 



a, ( + 2C - -Ai^ j 



9, ( C - ^Ai?) 



(3.37) 



-9H^^^ + 12H^C - 3C + QHCC - -H{AE){AE) 

3 

+2H{didjE){didjE) - {didjEf + ^(a^)^. 



(3.38) 



The difference of Eq. (3.38) and Eq.(3.37) gives precisely the torsion scalar (see Eq.(2.18)) 
up to second order in perturbation: 



T = 6i?2 + 12HC, - 12i?2$ - 4a-2(ai$) 



8.(c4ae) 



2a 



(3.39) 



+18^2^2 _ 2AH^C + 6C^ - 12HCC - 4:H {didj E) {didj E) 
+^H{AE){AE) - {d,d,Ef + i(A^)2. 

Note that Tq = 6H^ is simply the background torsion scalar for FLRW geometry with 
flat spatial section. Most calculations in the literature have the torsion scalar to be —6H^ 
instead. This difference by a minus sign can be directly checked to be the result of our choice 
of east coast sign convention (— , +, +, +) instead of west coast one (+, — , — , — ). 
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We can now expand the gravitational part of the action. Let us denote /o = /(Tq) and 
let prime denotes the derivative {df / dT)\-r=To ■ Then 



dtd'^x a"^ 



/o + (/o - 12//2/o)<^ + 3(/o - l2H^fo - ^Hf^ - Am^Hf'^X 
+QH\f'^ + 12^2 /o)^'' - 12^(/^ + 12^2 /o)^C + 3(/o - l2H^fo)K 



(3.40) 



where we have performed integration by parts and dropped the surface term. 

We must also expand the action for scalar field. We define the perturbation of scalar 
field as 



+ 



(3.41) 



where <j)Q is the background value of (j) and 6(j) is the perturbation of (j). Then, scalar field 
action is written up to the second order of perturbation as 



+ 



1 



^<^o-^+U'^o-^)(^ + 3C) + 



bl^ - V't 



l-V][m + 3C' - -{/\Ef ) + ($ + 3C)(<Ao 



(3.42) 



From the variation of the first order action with respect to ^ and 5(j), we can obtain 
the equations of motion of the background as 



/o - 12//2/0 



y = 0, 



(^0 + 3i?0o + V' = 



y = o, 



The first two equations give 



4F(/^ + 12//2/o)- 



(3.43) 

(3.44) 
(3.45) 

(3.46) 



We now proceed to the second order action. We shall switch back to the original variable 
ip from Q. The relation between these two variables is given by Eq.(3.23). Moreover, we use 
other choice of gauge where E = f3 = because the original gauge choice f3 = B = does 
not fix all gauge degrees of freedom. 

The second order action in the new gauge is easily obtained by replacing E with —aB 
in the corresponding equations with the original gauge choice'^ . Then, by integration by 



Gauge invariant combination including E is written as E — aB. See Eq.(3.14) and Eq.(3.16). 
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part, dropping the surface terms and using Eqs.(3.43)-(3.46), the second order action can be 
written as 



Si 



where 



2H 



3i72) 1 



H 



V'a + — v^n' 



lAB] 



-aAB + -r^a + 



MH + 3^2) 



-a 



-a 



+ 



+ 



IQf'^H^ 4/' 2Hf' 



a 



(3.47) 



a 



(3.48) 



We can see that there are no kinetic terms for ip and /3 and so we have only a single prop- 
agating degree of freedom. Moreover, a is of the same form as the gauge invariant variable 
in GR. As a result, the only difference from the corresponding second order perturbative 
expansion in GR appears as the replacement of 1/(2k) by an effective gravitational coupling 

3.3 Pseudoscalar Perturbation 

The cotetrads and metric for pseudoscalar perturbation are given by 



e"^(ix^ = dt, (3.49) 

e^^dx^ = 6aia [5ij + eijkdka] dx^ , (3.50) 

d d 

= ^aifl"^ + eijkdkd - 5ij{dkaf + {dia){dja)) (3.52) 

g^.udx^'dx'' = -dt^ + a2 {5,j + 6^j{^kaf - {dia){djd)) dx'dx^, (3.53) 
This leads to 

e = (1 + (didf) , (3.55) 

(3.56) 
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and the following components of torsion tensors: 

T\i = 0, (3.57) 

T% = 0, (3.58) 

T% = H5ij + e.jkdka + (du~a) (3^^) - {dfa) {did) , (3.59) 
jfc = eikididjo- - eijididkd + 5ik{dia){didja) 

-6ij{dia){didka) - {dkcr){didja) + {dja){didkd), (3.60) 

T^^^ = 2,H + 2{dia){d^d), (3.61) 

T% = {dfa){d,dfa) + {d,a){Aa). (3.62) 

We then compute the various contractions: 

T^,W = -9H^ - l2H{dia){did), (3.63) 

= -3H^ - m{dia){d^d) + a'^ [{didjdf - (Aa)^] . (3.64) 
Finally, we can obtain the torsion scalar 

T = 6H^ + m{did){did) + a'^ [{didjdf - (Aa)^] . (3.65) 
The second order action for gravitational field is 

Sgi = j dtcfx [fidiaf + SH f {dfa){did) + a'^/' [{8,8 faf - (Acr)^]] . (3.66) 



On the other hand, the second order scalar field action is 

1 



= / dtd?x 0? 



(3.67) 



Integration by parts followed by dropping the surface terms, and then using the back- 
ground equation Eq.(3.44), the summation of the above two actions S^f and gives zero. 
Thus, we see that there is no propagating pseudoscalar degree of freedom, up to second order 
in perturbation. 

3.4 Tensor Perturbation 

Moving on to tensor perturbation, the cotetrads and metrics are given by 



e^^dx^' = dt, (3.68) 

e^dx" = 6ai (<5,j- + ) dx^ , (3.69) 

8 8 

^a'^-Q^ = a-^^aj{Sij - hij + hikhkj)-^, (3.71) 

gfj_i,dx'^dx'^ = —dt^ + a^{6ij + 2hij + hikhkj)dx'^dx\ (3.72) 
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We thus have 



e = a'^ ( 1 - -hij 



The components of the torsion tensor are: 



rt 
i 



0, 



T* • = 



r 



tj ■ 

n 
jk 



H6ij + hij 

■ djhik — dkhij — hiidjhik + hudkhi 



— 3i? hij hij, 
T^i^ = —hjkdihjk + hjkdjhki- 

We proceed to calculate in the same way we did before 

= -9H^ + 6Hh,jkj, 

fiu-'- 13^9 pay 9 ^2 

= -3H^ - ill + '^Hhijhij + [{dihjkf - {dihjk){djhik)\ 
and obtain the torsion scalar 

T = m^- h% - AHh.hij + [{d.hjkf - {dihjk){djhik)\ . 
Then, the second order action for tensor mode is found to be 



(3.74) 
(3.75) 



(3.76) 
(3.77) 
(3.78) 
(3.79) 
(3.80) 
(3.81) 



(3.82) 
(3.83) 

(3.84) 



S. 



9,2- 



-f'^hlj - AHf^hijk, + a-2/o [{dih.k? - {^^hjk){^,h,k)\ - -kh% 



On the other hand, the second order contribution from the scalar field action is 

dt(fx 



(3.85) 



(3.86) 



We can simplify the summation of the above two actions by performing integration by parts, 
dropping the surface terms and using background equation (3.44). 
We have 



5j 



dtdPx a^-f'o) hfj - a-\d^hjk) 



(3.87) 



As with the scalar perturbation, the only difference of the tensor perturbation from that of 
its counterpart in GR is the overall factor. While in GR the overall factor is l/(2«;), in f{T) 
gravity we have the time-dependent coefficient — /q- 
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3.5 Vector and Pseudovector Perturbations 

Recall that the pseudovector degree of freedom is coupled to the vector degrees of freedom. 
Nevertheless by gauge freedom, we can fix the gauge as 1^ = 0, so that we need not worry 
about the pseudovector degrees of freedom. 

The cotetrad and the metric are given by 



e ^dx^ = dt + auidx^ 



e°-^dx^ = 6ai [vidt + a {6ij + diWj + djWi) dx^] , 



,. d 



d d 

(1 + UiVi) — + a^-^[-Vi + Vj{diWj + djWi)]-^, 

d 

{-Ui + Uj{diWj + djWi)} — 



d 

+a"^ {5ij - {diWj + djWi) + {diWk + dkWi){djWk + dkWj) + UiVj} 



(3.88) 
(3.89) 

(3.90) 



(3.91) 



Qp^pdx^dx'^ = — (1 — vf)dt^ + 2a[—Ui + Vi + Vj{diWj + djWi)]dtdx^ 

[6ij + 2{diWj + djWi) - UiUj + {diWk + dkWi){djWk + dkWj)] dx^dx^ , (3.92) 



_d d_ 

dx>^ dx^ 



-(1 + 2uiVi - uf ) 



dt 



+2a ^[-Ui + Vi + 2uj{diWj + djWi) - Vj{diWj + djWi)] 
+a~'^[6ij - 2{diWj + djWi) + 2uiVj - ViVj 



d__d_ 

dt dx^ 



+3{diWk + dkWi){djWk + dkWj)] 



d d 

dx^ dx^ 



Then, we have 



(3.93) 



1 - UiVi - -{diWj + djWi 



(3.94) 
(3.95) 



From the components of the torsion tensor: 



T' 



(3.96) 
(3.97) 



(3.98) 



a[ui - Uj{diWj + djWi)] + ujdiVj, 
a[{diUj - djUi) - UkdkidiWj - djWi)], 
T\j = H6ij + {diWj + djWi) - a~^djVi- Viiij 

-{diWk + dkWi){djWk + dkWj) + a~^{diWk + dkWi)djVk, 
Tjk = di{djWk - dkWj) - Vi{djUk - dkUj) - {diWi + diWi)di{djWk - dkWj), (3.99) 
r'^^ = 3H - ViUi - {diWj + djWi){diWj + djWi) + a~^{diWj + djWi)diVj, (3.100) 
-aiii — Awi + auj{diWj + djWi) — ujdiVj — Vj{diUj — djUi) 
-{dkWj + djWk)dk{diWj - djWi), (3.101) 
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we can compute 

Tf.T^g'"' = -9H'^ + + muiiii + 9H'^uj - im'^mvi + 2a~^UiAwi 

+6Ha~^UiAwi — 6Ha~^diVj{diWj + djWi) — 6Ha~^ViAwi 
+6H{diWj + djWi){diWj + djWi) + a~2(A^^;^)^ (3.102) 



1 

4 

= ~2^^ + ~^ Huiiii + -j^H'^u^i + -^a~'^{diUj - djUif - SH'^UiVi 
-^a~'^{diVj)'^ + Ha^^UiAwi + a^^{diVj){diWj + djWi) 
—Ha~^{diVj){diWj + djWi) — Ha~^ViAwi — -{diWj + djWif' 
+H{diWj + djWi){diWj + djWi) + ^a~'^[di{dkWj + djWk)]'^, (3.103) 

3 1 3 

-^a~'^{diVj){djVi) + Ha'^UiAwi + a~^{diVj){diWj + (9jmj) 

— H {diVj){diWj + djWi) — Ha~^ViAwi — -{diWj + djWi)^ 
+H{diWj + djWi){diWj + 5jWi) 

+ + djWk)][dj{dkWi + SjWfc)]. (3.104) 

Consequently, the torsion scalar is given by 
T = 6H^ — AHuiiii — QU'^uf — -a~'^{diUj — djUi)^ + a~'^diVj{diUj — djUi) 

+12H'^UiVi — -o.~'^ (diVj)'^ — -a~'^ {diVj){djVi) — 4Ha~^UiAwi — 2a^^UiAwi 

+2a~^ {diVj){diWj + djWi) + 4Ha~^ {diVj){diWj + djWi) + 4Ha~^ViAwi 
-{diiiij + djWi)'^ - AH{diWj + djWi){diWj + djWi) - a~'^{Awif 

+ ^a~'^[di{dkWj + djWk)f + ]^a~'^[di{dkWj + djWk)][dj{dkWi + diWk)]. (3.105) 
Therefore, the second order action for vector mode is given by 

SI2 = j dtd^x [U-um + w^Awi) + f(^62T] . (3.106) 

where is the second order part of T. The second order part of the scalar field action is 

1 



S^2= I dtd-'x a' 



^(pliuiVi -uf + WiAwi) + V{uiVi - WiAwi) 



(3.107) 



We can again simplify the above actions via integrating by parts, dropping the surface terms 
and using background equations Eqs.(3.43)-(3.46), we can finally obtain the second order 
total action as 

S2 = j dtd^x -af'^ivi — Ui — 2awi)A{vi — Ui — 2awi). (3.108) 
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From this we can see that there is no propagating degree freedom of vector mode up to second 
order in linear perturbation. 

4 Discussion 

In conclusion, we examined second order perturbation of f{T) theory of gravity and found 
that there is no scalar, pseudoscalar, pseudovector or vector modes of gravitational propaga- 
tion. We also found that the scalar and tensor perturbations have exactly the same form as 
their counterparts in usual general relativity with scalar field, except that the factor reduced 
Planck mass squared Mpj that occurs in the latter has now been replaced by an effective 
time-dependent gravitational coupling —2{df / dT)\T=To : where Tq is the background torsion 
scalar. 

We nevertheless emphasize that our current calculation only reveals dynamics of the 
propagating modes in f{T) gravity in a highly symmetric background of FLRW geometry, 
which is although useful for the purpose of studying some aspects of cosmological properties, 
does not reveal the full features of what f{T) theory has to offer. In general /(T) gravity has 
three extra degrees of freedom in nonlinear analysis [11], and the extra degrees of freedom 
might be excited in spacetime with less symmetry. 

The situation is very similar to that in nonlinear massive gravity theory [25] in which, 
while nonlinear analysis shows that there are generically five gravitational degrees of freedom 
[26], in the second order action on open FLRW background we can catch only two tensor 
degrees of freedom [27] [28]. Since the hidden degrees of freedom cause nonlinear instability 
in the case of nonlinear massive gravity [29], we believe that one cannot conclusively and 
confidently decide on the issues of stability in /(T) theory as long as one does not know 
the behaviors of all the degrees of freedom that might manifest themselves in more generic 
conditions, and this is an important issue in cosmology since the real universe is certainly 
not completely homogeneous and isotropic. 

A fully nonlinear analysis for f{T) gravity is of course difficult to carry out. It is thus 
an interesting problem whether we can see the effects of these degrees of freedom at linear 
perturbation level of anisotropic model of cosmology [30], e.g. Bianchi Type-I geometry 
[31], since such perturbation can be viewed effectively as nonlinear perturbation on FLRW 
cosmology. See e.g. [29] for such analysis in the context of nonlinear massive gravity. We 
will address these interesting issues elsewhere. 
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